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In the cosmological paradigm, cold dark matter (DM) dominates the mass content of the Universe and is present 
at every scale. Candidates for DM include many extensions of the standard model, with a weakly interacting 
massive particle (WIMP) in the mass range from ^^10 GeV to greater than 10 TeV. The self-annihilation or 
decay of WIMPs in astrophysical regions of high DM density can produce secondary particles including very 
high energy (VHE) gamma rays with energy up to the DM particle mass. VERITAS, an array of atmospheric 
Cherenkov telescopes, sensitive to VHE gamma rays in the 85 GeV - 30 TeV energy range, has been utilized 
for DM searches. The possible astrophysical objects considered to be candidates for indirect DM detection are 
VERITAS dwarf spheroidal galaxies (dSphs) of the Local Group and the Galactic Center, among others. This 
presentation reports on our extensive observations of these targets and constraints of the dark matter physics 
from these objects, including the methodology and preliminary results of a combined DM search of five dSphs. 


1. INTRODUCTION 

The search for Standard Model (SM) particles re¬ 
sulting from the annihilation of Dark Matter parti¬ 
cles provides an important complement to that of di¬ 
rect searches for DM interactions and accelerator pro¬ 
duction experiments. Among the theoretical candi¬ 
dates for the DM particle [H , weakly interacting mas¬ 
sive particles are well motivated since they naturally 
provide the measured present day cold DM density 

. Candidates for WIMP dark matter are present in 
many extensions of the SM of particle physics, such 
as supersymmetry (SUSY) Q or theories with extra 
dimensions Q. In such models, the WIMPs either 
decay or self-annihilate into standard model particles, 
most of which produce either a continuum of gamma 
rays with energies up to the DM particle mass, or 
mono-energetic gamma-ray lines. 

Attractive targets for indirect DM searches are 
nearby massive objects with high inferred DM density 
which are not expected to be sources of VHE gamma 
rays. The Galactic Center is likely the brightest source 
of gamma rays resulting from DM annihilations, how¬ 
ever the detected VHE gamma-ray emission is coinci¬ 
dent with the supermassive black hole Sgr A* and a 
nearby pulsar wind nebula @ , motivating searches for 
DM annihilation in the Galactic Center halo where the 
VHE gamma-ray background is expected to be signifi¬ 
cantly lower @ . Galaxy Glusters have a large DM con¬ 
tent. However, they are extended for VERITAS, and 
the possibility exists of a VHE background from con¬ 
ventional processes a a, although not yet detected. 
Galactic DM sub-halos would appear as unidentified 
objects (UNIDs) without multi-wavelength counter¬ 
parts in Fermi-LAT data. If a Fermi UNID were 
detected in VHE, it could potentially be from DM. 
Dwarf spheroidal galaxies (dSphs) are additional at¬ 
tractive targets for DM searches. Dwarf spheroidal 
galaxies are relatively close (^50 kpc), and have a 
low rate of active or recent star formation, which sug¬ 


gests a low background from conventional astrophys¬ 
ical VHE processes 

The following sections describe the status of obser¬ 
vations and data analysis of each of the DM targets 
described above as of fall 2014, followed by conclu¬ 
sions and a discussion of the future of the VERITAS 
DM program. 


2. SUB-HALOS 

Recent cosmological N-body, high-resolution simu¬ 
lations M indicate that DM halos are populated with 
a wealth of substructures m- Because of tidal dis¬ 
ruption near the Galactic disk, most of the sub-halos 
are thought to survive at high Galactic latitude. The 
lack of material in these regions prevents the DM over¬ 
densities from attracting enough baryonic matter and 
trigger star formation. DM clumps would therefore 
be invisible to most astronomical observations from 
radio to X-rays. DM structures residing in the the 
Milky Way halo can be nearby the Sun and there¬ 
fore have a bright gamma-ray annihilation signal [l^ . 
These clumps would likely be only visible at gamma- 
ray energies and therefore may not have shown up in 
astronomical catalogs yet. Since gamma-ray emission 
from DM annihilation is expected to be constant, DM 
clumps could then appear in all-sky monitoring pro¬ 
grams [l^ done at gamma-ray energies. These can 
be best provided by the Fermi-LAT instrument. Very 
likely, the distinct spectral cut-off at the DM particle 
mass is located at energies too high to be measur¬ 
able by Fermi within a reasonable timescale (see, e.g. 
the WIMP mass lower limits in El) and can only be 
detected by ground-based telescopes, such as VERI¬ 
TAS. Furthermore, detection of this spectral cut-off 
at the same energy in multiple objects would stand 
as a visible signature of DM. The Second Fermi-LAT 
Catalog (2FGL) contains 1873 high energy gamma- 
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Table I Preliminary DM Sub-halo Results. Flux upper limits are given in units of Crab Nebula flux. 


2FGL Name 

Exposure (hrs) 

Significance (a) 

Excess Counts 

Ft, (GeV) 


J0312.8-k2013 

9.7 

-1.5 

-26 ± 17 

220 

< 0.9% 

J0746.0-0222 

9.1 

-0.9 

-15 ± 16 

320 

< 1.1% 


ray sources detected by the LAT instrument after the 
first 24 months of observations. For each source, posi¬ 
tional and spectral information are provided as well 
as identification or possible associations with cata¬ 
loged sources at other wavelengths. Although Fermi- 
LAT has a good angular resolution, a firm identifica¬ 
tion based on positional coincidence alone is not al¬ 
ways feasible. Thus, 576 sources in the 2FGL lack 
any clear association. These are the so-called unas¬ 
sociated Fermi objects (UFOs), a population among 
which DM clumps might be represented [l^ . In or¬ 
der to extract possible DM clump candidates out of 
the 2FGL UNIDs, we adapt the selection criteria from 
selecting sources by requesting them: 

• to lie outside the Galactic Plane, 

• to be non-variable, 

• to exhibit a power law spectra, and 

• to not have possible counterparts. 

The original list obtained from the 2FGL catalog 
is then filtered to select only sources observable with 
VERITAS with a maximum zenith angle at culmi¬ 
nation of 40°, in order to pursue the lowest energy 
threshold. Additionally, an estimate of required ob¬ 
servation time for a 5cr detection, dubbed detectability, 
is computed based on a 2FGL Gatalog flux extrapo¬ 
lation to the VERITAS energy range. 

The preliminary results of the VERITAS observa¬ 
tions shown in Table 1 are in tension with the ex¬ 
trapolation of the gamma-ray spectra from the Fermi- 
LAT to very high energies. Additional data from these 
UNIDs and others by VERITAS and other Cherenkov 
telescopes may completely completely rule out a di¬ 
rect extrapolation of the Fermi-LAT which would give 
strong DM model constraints or potentially detect a 
DM signature, provided that they are truly without 
counterparts at other wavelengths. 

3. GALAXY CLUSTERS 

Clusters of galaxies are the largest viralized objects 
in the Universe, with typical sizes of a few Mpc and 
masses on the order of 10^"* to 10^® M 0 . Most of 
the mass 80%) is dark matter, as indicated by 
galaxy dynamics and gravitational leasing [T^. Aside 
from DM annihilation, it is possible to have gamma- 
ray emission from cosmic-ray interactions, producing 



Figure 1: Dark Matter velocity-averaged cross-section 
limits from the Coma galaxy cluster. Figure taken from 

0 . 


neutral pions Q , or inverse Compton of ambient pho¬ 
tons 0 . 

VERITAS has taken 18.6 hours of dedicated ob¬ 
servations of the Coma cluster between March and 
May 2008. The Coma cluster is a close (z=0.023) 
and massive (M~ cluster which has been 

thoroughly studied across all wavelengths. The stan¬ 
dard analysis of the Coma cluster data using point- 
source cuts for the core of the cluster yielded 17 ex¬ 
cess counts, with a significance of 0.84cr, indicating a 
non-detection. Upper limits of 0.83% of the Crab Neb¬ 
ula flux were placed for the core of the Coma cluster 
with 95% confidence, assuming a powerlaw spectral 
index of -2.3. With the absence of a signal from the 
Coma cluster, limits of the velocity-averaged cross- 
section for DM annihilation were placed at 0(10“^^), 
as shown in Figure 1 M- 

An archival VERITAS galaxy cluster search is also 
in the works, looking for galaxy clusters that have hap¬ 
pened to overlap in the same FOV as other targeted 
observations. ROSAT and SDSS galaxy cluster cat¬ 
alogs are being used to cross-check with other VER¬ 
ITAS observations. Most notably M87 in the Virgo 
cluster of g alaxies and NGC 1275 in the Perseus 
cluster | 2 l| . 
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Figure 2: Skymap of the galactic center region using a 
subset of the VERITAS data. DM signal and background 
regions are indicated north and south of SgrA*. 


4. GALACTIC CENTER 


The center of our galaxy, SgrA*, is a strong VHE 
source, along with several other VHE sources nearby 
and possible diffuse emission (2^, making Dark 
Matter detection in that region a complicated, but 
not impossible, prospect. 

The Galactic center was observed by VERITAS 
in 2010-2014 for ^ SOhrs (good quality data) at 
zenith angles of z = 60 — 66 deg (average threshold 
of Ethr — 2.5 TeV). The higher effective areas due to 
the large zenith angle observations make the VERI¬ 
TAS observations now the most sensitive instrument 
for the Galactic Center region above 2 TeV. The detec¬ 
tion of SgrA* by VERITAS and VHE emission in the 
region through conventional processes are discussed in 
greater detail elsewhere in these proceedings [^ . 

The DM signal and background regions for the 
Galactic center region will use arc-shaped regions 
north and south of the Galactic plane to avoid dif¬ 
fuse emission and VHE sources, as shown in Figure 
2. The VERITAS observations were accompanied by 
off-source observations of a field located in the vicin¬ 
ity of the Galactic center region (with similar zenith 
angles and sky brightness) without a known TeV 7 - 
ray source. These observations are used to study the 
background acceptance throughout the field of view 
and will assist in the identification of a diffuse 7 -ray 
component surrounding the position of the Galactic 
center. 

The DM search for the Galactic center region is still 
in the preliminary stages. Work is currently underway 
for computing J factors for the signal and background 
regions for the arc-shaped regions described above. 


5. DWARF GALAXIES 


Dwarf spheroidal galaxies (dSphs) best meet the cri¬ 
teria for a clear and unambiguous detection of dark 
matter. They are gravitationally-bound objects and 
contain up to 0(10^) times more mass in DM than 
in visible matter [1]. As opposed to the Galactic cen¬ 
ter, globular clusters and clusters of galaxies, dSphs 
present the clear advantage of being free of any sig¬ 
nificant astrophysical emission. Their high Galactic 
latitude and relative proximity to Earth (^ 50 kpc) 
make them very good targets for high signal-to-noise 
detection. 

Between the start of full VERITAS array operations 
and Summer 2013, five dSphs have been observed with 
VERITAS: Segue 1, Ursa Minor (UMi), Draco, Bootes 
1, and Willman 1. The VERITAS collaboration has 

g viously published a 48 hour exposure on Segue 1 
and ^^15 hour exposures on the other four men¬ 
tioned here [23 ■ Deeper exposures on Segue, UMi and 
Draco have been taken after these publications. To 
obtain the lowest possible energy threshold for DM 
searches, looser cuts optimized apirori on soft spec¬ 
tral VHE sources were used for the collective data set. 
The combination of looser cuts and deeper exposures 
revealed certain systematic effects in the cosmic-ray 
(GR) background data. The first is a gradient in the 
VERITAS cameras dependent on the zenith angle of 
observations across the FOV. The second systematic 
effect results from bright stars in the VERITAS FOV 
that cause the high voltage to pixels in the cameras to 
be suppressed. Both of these systematic effects have 
been corrected for and the results are summarized in 
Table 2. The ‘crescent’ background method which also 
developed for dSph analysis [27| was also used for the 
Table 2 results. 

The first background systematic effect, relating to 
the zenith gradient, was corrected using a zenith- 
dependent acceptance map. The standard VERITAS 
analysis uses only a radially-dependent acceptance, 
i.e. the angle between the reconstructed event direc¬ 
tion and the array pointing direction. Measuring the 
gradient utilized a map that is the ratio of the num¬ 
ber of all reconstructed events in a sky map within 
a given search radius (defined as 0.17 degree in this 
work) to the (radial only) acceptance in that same bin, 
a parameter we will refer to as flatness in the rest of 
this work. If the acceptance adequately describes the 
GR background, then excluding any stars or known 
VHE sources, it should not correlate with any exter¬ 
nal parameters. However, a strong correlation was 
seen with the mean zenith angle of each reconstructed 
event position in the skymap bin. This correlation in 
the skymap bins is shown in Figure 3. This gradi¬ 
ent is corrected in the data by fitting the correlation 
with a fourth-degree polynomial and using that to re¬ 
weight the acceptance map. The a parameter from 
Li & Ma equation 17 is then re-calculated [ 2 ^, [ 2 ^. 
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Table II Preliminary DM DSph Results 


DSph Name 

Exposure (hrs) 

Significance (cr) 

Excess Connts 

Etr (GeV) 

F^^^cn^C.U.) 

Segue 1 

92.0 

0.7 

94.4 ± 134.1 

150 

< 0.4% 

Ursa Minor 

59.7 

-0.1 

-7.2 ± 68.5 

290 

< 0.3% 

Draco 

49.9 

-1.0 

-73.2 ± 69.1 

220 

< 0.3% 

Bootes 1 

14.0 

-1.0 

-38.5 ± 36.7 

170 

< 0.3% 

Willman 1 

13.7 

-0.6 

-28.7 ± 46.2 

180 

< 1.0% 



A Zn (deg) 


Figure 3: Scatter plot of the flatness parameter on the 
y-axis, and the mean zenith difference between the array 
tracking direction and event reconstruction direction on 
the x-axis for the Segue 1 data summarized in Table 2. A 
fit of this scatter to a fourth-degree polynomial is shown 
in red. 
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Figure 4: Example of the HFit and standard Hillas event 
characterization. The HFit shower image is the blue 
ellipse outline, while the standard Hillas moment analysis 
is the red ellipse outline. Figure from [^ . 


It should be noted that the difference of the adjusted 
value of a to the non-adjusted value is typically less 
than 1%. However, the difference to the y-ray excess 
and significance becomes larger over time as statistics 
accumulate. 

The second background systematic effect is due to 
“holes” that are seen in the data relating to bright 
stars in the FOV that would trip the high voltage 
of camera pixels or raise cleaning thresholds due to 
higher night-sky background levels. Missing pixels 
would in turn effect both energy and gamma-ray po¬ 
sition reconstruction. A method of using a 2D Gaus¬ 
sian likelihood fit to each shower image is utilized here, 
called HFit js^]. Standard VERITAS analysis uses the 
moments of the shower images, commonly referred to 
as Hillas parameters. By using the 2D elliptical fitting 
to each image, missing data from disabled or broken 
PMTs are effectively interpolated around, as are im¬ 
ages truncated by the edge of the cameras, as shown 


in Figure 4. This has been shown to reduce both the 
size and the depth of the holes due to bright stars 
seen in the data, including but not limited to the B 
magnitude 3.4 star Eta Leonis which is located 0.68° 
away from the center of the Segue 1 dSph. The ef¬ 
fectiveness of HFit on a independent data sample is 
shown in Figure 5, which shows the apparent surface 
brightness in the CR background (which is in reality a 
deficit for reasons described above) at a star location 
in the FOV of the blazar RGB JI058-I-564 (Merak, 2.4 
B magnitude). 

Work is currently underway to utilize the data for 
the previously published papers plus additional data 
for a combined DM physics result. This result will use 
the methodology developed by Geringer-Sameth et al. 
(sTjl to utilize both the individual energy and event 
reconstruction information as well as astrophysical “J 
factors” from a generalized NFW profile by Geringer- 
Sameth and Walker [s^- 
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Theta square plot (Wobble) 
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Figure 5: Apparent surface brightness of background 
cosmic ray events as a function of angle 9 from the bright 
star, Merak, in the FOV of RGB J1058+564 using both 
the HFit algorithm and the standard Hillas shower image 
characterization. The VERITAS data analysis typically 
excludes regions around bright stars or known VHE 
sources from cosmic-ray background characterization. 

The default radius for background exclusion region for 
Merak is also shown as the dashed blue line. 

6. CONCLUSIONS 

New DM publications from the VERITAS collabo¬ 
ration are forthcoming: the combined analysis of the 
dSphs should be publically available within the next 
six months, which promises to be the most robust re¬ 
sult of any DM result in VHE gamma rays so far, while 
DM results for the Galactic Center, Eermi UNIDs and 
the archival galaxy cluster search should be ready on 
longer timescales. New analysis techniques are be¬ 
ing developed by the VERITAS collaboration which 
promise large gains to our DM sensitivity, as an exam¬ 
ple an extended analysis of the dSphs which would in¬ 
corporate longer tails of the DM density profile, which 
would in turn give a boost to the J factors by as much 
as a factor of ~2. A combined analysis with Eermi- 
LAT, or other 7 -ray instruments could potentially pro¬ 
vide a boost to DM sensitivity. 

The VERITAS collaboration has a historical com¬ 
mitment to substantial DM observations and plans to 
do so in the foreseeable future. Recently, a new long¬ 
term plan for VERITAS has gone into effect, which 
has a significant fraction (15-20%) of the total dark 
observation time dedicated to DM observations. The 
focus of this long-term plan is dSphs; however the 
galactic center, Eermi UNIDs and galaxy clusters will 
not be completely ignored. If executed consistently 
over the expected lifetime of VERITAS, these obser¬ 
vations will form the basis of an important and unique 
contribution to the field of indirect DM detection. 
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